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Abstract. First, we recall the algebro-geometric method of con- 
struction of finite field valued solutions 

Z 3 3 (ni, 712,713) 1 — ► r(ni, 712,713) G F 9 , 

of the discrete KP equation (T — shift operator in 7ij variable) 

(Tit) • (T 2 T 3 r) - (T 2 t) • (T 3 T x t) + (T 3 r) • (TlT 2 t) - 0, 

and next we perform a reduction of dKP to the discrete ID Toda 
equation 

(T^Tst) • (T^-V) - (TfV) • (Tr) + r • r = 0. 

This gives a method of construction of solutions of the discrete ID 
Toda equation taking values in the finite field ¥ q . 



1. Introduction 

Discrete integrable systems focus much attention for at least two 
reasons. They inherit developments of the theory of continuous inte- 
grable models on the one hand, and give rise to new ideas on the other. 
Recently this "twofold" nature was exploited for the construction of 
integrable cellular automata. 

The approach is based on transferring the standard algebro-geometric 
method of construction of solutions from the complex field C to a 
finite field case. In spite of apparent differences between complex and 
finite fields (see [TT]) the theory of algebraic functions over finite field 
El provides a powerful framework and makes this transfer possi- 
ble. As a consequence of algebraic approach it appears that relevant 
theorems (as the Riemann-Roch theorem) and formulas are " quite the 
same" as in the complex field case but have a different meaning. For 
instance, the concept of a continuum limits (in expansions etc.) is not 
clear as far as finite fields are regarded. 
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For a given difference equation, its finite field valued solution de- 
fines a cellular automaton. So far there have been constructed mul- 
tisoliton solutions for the discrete 2D Toda lattice (the Hirota's bilin- 
ear difference equation) jH], the discrete Kadomtsev-Petviashvili (dKP) 
and discrete Korteweg-de Vries (dKdV) equations [3| and also algebro- 
geometric solutions for dKP jl] and dKdV |2j. In this article we extend 
previous results for the discrete ID Toda equation (and ipso facto for 
a certain discrete version of the sine-Gordon equation, see remark at 
the end of section 13. 1|) . 

The paper consists of three main parts. In Section |21 we recall the 
finite field version of the algebro-geometric construction of solutions 
of the discrete KP equation. In Section |3] the algebro-geometric re- 
duction to the discrete 1 dimensional Toda equation is performed and 
illustrated by examples in Section 

2. An algebro-geometric approach to dKP equation 

2.1. General construction for the dKP equation. For the gen- 
eral construction |3j we need an algebraic projective curve C/K (or 
simply C), absolutely irreducible, nonsingular, of genus g, defined over 
the finite field K = ¥ q with q elements. C(K) denotes the set of Ir- 
rational points of the curve. K denotes the algebraic closure of K, i.e., 
K = IJ^Fg*, an d C(K) denotes the corresponding infinite set of Ir- 
rational points of the curve. Denote by Div(C) the abelian group of 
the divisors on the curve C. The action of the Galois group G(K/K) 
(of automorphisms of K which are identity on K) extends naturally 
to action on C(K) and Div(C). A field of K-rational functions on 
the curve C is denoted by K(C). A vector space L(D) is defined as 
{/ G K(C) | (/) > -D}, where D G Div(C) and (/) = £ PeC ord P (f)-P 
is the divisor of the function / G K(C). 
On the curve C we choose: 

(1) four points A , A t G C(K), % = 1, 2, 3, 

(2) effective K-rational divisor of order g, i.e., g points _E> 7 G C(K), 
7 = 1, . . . ,g, which satisfies the following K-rationality condi- 
tion 

Vcr G G(K/K), a(fi 7 )= J B y , 

(3) effective K-rational divisor of order N, i.e. N points C a G C(K), 
a = 1, . . . , N, that satisfies the K-rationality conditions 



Va G C7(K/K), a(C a ) = C a ,, 
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(4) N pairs of points Dp,E p e C(K), (3 = 1,...,N, that satisfy 
K-rationality conditions 

W e G(K/K) : a{{Dp,E p }) = {D p ,,E p ,}. 

We assume that all the points used are distinct and in general posi- 
tion. In particular, the divisors ^ 7 =i -^t> Ef=i^«) an d D(ni, ri2, n 3 ), 
defined below, are non-special. 

Definition. Fix K-rational local parameter t at A . For any integers 
ni,n2,n 3 e Z let divisor D(ni, 712,713) be of the form 

3 g N 

D(n u n 2 , n 3 ) = n^A, - Aj) + S 7 + C a . 

i=l 7=1 ct=l 

The function ^(rii, ri2, 713) (called a wave function) is a rational function 
on the curve C with the following properties 

(1) the divisor of the function satisfies (ip) > —D, i.e. ip G L(D), 

(2) the function ip satisfies N constraints 

i/j(n 1 ,n 2 ,n 3 )(D f3 ) = ^(ni, n 2 , "-3) (#/?), /3 = 1, ... ,7V. 

(3) the first nontrivial coefficient of its expansion in to & t Aq is 
normalised to one. 

Existence and uniqueness of the function ^(ni, 712, n 3 ) is due to appli- 
cation of the Riemann-Roch theorem with general position assumption 
and due to normalisation. Moreover, the function ip(ni, n 2 , n 3 ) is DC- 
rational, which follows from K-rationality conditions for sets of points 
in their definition. 

The next step of the construction is to obtain linear equations for 
wave functions. The full form of such equation is in case when the 
pole of ip(ni,n 2 ,n 3 ) at A is of the order exactly (rii + n 2 + n 3 ) and 
respective zeros at A { are of the order rij, for i = 1,2, 3. We will 
call this case generic. Having fixed IC-rational local parameters t« at 
Ai, % = 1,2,3, denote by Ck\ n i-> n 2, ^3), i — 1,2,3, the K-rational 
coefficients of expansion of ip(ni, n 2 , n 3 ) at A iy respectively, i.e., 

00 

V>(ni, n 2 , n 3 ) = C^K, n 2 , n 3 )tj, i = 1, 2, 3. 

fe=0 

Denote by Tj the operator of translation in the variable % — 1,2, 3, 
for example, T 2 ip(ni, n 2 , n 3 ) = ^(ni, ^2+1, ^3)- The full linear equation 
is of the form 

(1) T ^-T^ + ^^ = 0, i^j, i,j = 1,2,3. 

Co 
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It follows from observation that Tiip — Tjip e L(D), hence it must be 
proportional to wave function ip. Coefficients of proportionality can be 
obtained from comparison (the lowest degree terms) of expansions of 
left and right sides of (JJJ at the point A4. 

Remark. When the genericity assumption fails then the linear prob- 
lem (£Q) degenerates to the form Tiip = ip or even to = 0. 

Notice that equation (JTJ) gives 



(2) ^ = --%-> **i< «' = 1 . 2 . 3 - 

so so 

Define 

Pi = {-ljZi« n '$>, 1 = 1,2,3, 

then equation (J2J implies existence of a K-valued potential (the r- 
function) defined (up to a multiplicative constant) by formulas 

(3) — = Pi, i = l,2,3. 

r 

Finally, by " cyclic" use of equations JTJ) and (J2J) one can get the condi- 
tion 

(4) T,p 1 _np 1 + np l = ^ 

Pl Pl P2 

which, written in terms of the r-function, gives the discrete KP equa- 
tion [Hj called also the Hirota equation 

(5) (T lT ) (T 2 r 3 r) - (T 2 r) (T 3 T lT ) + (T 3 r) (T^r) = 0. 

Remark. Equation (j3J) can be obtained also from expansion of equation 
(JTJ at Ak, where k = 1, 2, 3, k ^ i, j. 

Absence of a term in the linear problem (0), due to the Remark 
above, reflects in absence of the corresponding term in equation (jSJ). 
This implies that in the non-generic case, when we have not defined 
the r-function yet, we are forced to put it to zero. 

2.2. Explicit formulas for projective line. Denote by t standard 
parameter on projective line P(K) and choose A = oo. Then we have 
the explicit formulas: vacuum wave functions 

^° = (t-A 1 ) n i(t-A 2 ) n °(t-A 3 ) n * 

and vacuum r-functions 

r° = (A 1 - A 2 ) nin2 {A 1 - A 3 ) nins (A 2 - A 3 ) n2n \ 
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Non vacuum r-function can be constructed by (see [Hj) 
(6) T = T°det(j) A (D,E), 

where we denote by <p A (D, E) the N x N matrix whose element in row 
j3 and column a is 

[0° A (A E)] ap = fijpp) -C(Ep), a, (3 = 1, . . . , N. 

and auxiliary vacuum wave functions 0° , a — 1, . . . , N have the form 
, _ 1 {t-A 1 ) n -{t-A 2 )^{t-A^ 



a 



t-C a (C a -A 1 )<*{C a -A 2 ) n *{C a -A 3 ) 
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Remark. Let L = ¥ q i be a field of rationality of all the points used in 
the construction. Then r-function is periodic in n\, n 2 and n 3 with 
periods being divisors of q l — 1, which is the order of the multiplicative 
group L*. 

Algebro-geometric solutions can be obtained using the Jacobian of 
the curve C. Construction of KP cellular automata from hyperelliptic 
curves was presented in j3] and of KdV cellular automata in [2] . 

3. Reduction to the discrete ID Toda equation 

3.1. Reduction to the discrete ID Toda equation. The discrete 
ID Toda equation 

(7) (T^-V) (Tf ^r) - (T lT ) (Tf V) + r r = 0. 
can be obtained from the dKP equation 

(T ir ) (T 2 T 3 r) - (T 2 r) (T 3 T ir ) + (T 3 r) (TxT 2 r) = 0. 
by imposing the constraint 

(8) TiT 2 T 3 -V = jt, 

where 7 is a non- zero constant. The reduction (jHJ) can be expressed at 
the level of the wave function ip as follows. 

Lemma. Assume that on the algebraic curve C there exists a rational 
function h with the following properties 

(1) the divisor of the function is (h) = —A + A x + A 2 — A 3 , 

(2) it satisfies N constraints: h(Dp) = h(Efs), (3 — 1, . . . , N , 

(3) the first nontrivial coefficient of its expansion in the parameter 
to at Aq is normalised to one. 

Then the wave function ip satisfies the following condition 

riT 2 T 3 -V = H. 
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Proof. Orders of poles and zeros of functions from both sides of equa- 
tion (J3.1j) are the same, so these functions are proportional. Normal- 
isation conditions (item |3] in Lemma and item |3] in Definition) imply 
they are equal. □ 

Remark. Equation (JJJ) is identical with a discrete version of sine-Gordon 
equation considered in [TUj . It is different from the standard form 
obtained by Hirota. As a consequence, cellular automata proposed 
here are different from those investigated in j3]. 

3.2. A r-function for the discrete ID Toda equation. Using re- 
duction condition for the wave function ip we can derive r-function 
satisfying the discrete ID Toda equation. 

Proposition. Let h be the function as in Lemma \3. 11 Let 8%, 5 2 and 

5 3 denote the respective first coefficients of local expansions of h in 
parameters t±, t 2 o,nd t 3 at A\, A 2 and A 3 , 

h = t 1 (5 1 + ...), h = t 2 (6 2 + ...), h= ^(6 3 + ...)• 

£3 

Then the function 

(g) f _ T fi- n i( n i- 1 )/ 2 ^_fi^-n 2 (n 2 -l)/2fin 3 (n 3 ~l)/2 

satisfies the discrete ID Toda equation (|7J). 

Proof. Expansions of equation (|3.1|) at points A%, A 2 and A 3 give 
T\T 2 T 3 x p x = Sxpt, T 1 T 2 T 3 1 p 2 = -5 2 p 2 , T 1 T 2 T 3 ~ 1 p 3 = 5 3 p 3 . 
The functions 

Pi = 5^ Pl , p 2 = (-5 2 )-> 2 , p 3 = 5 3 n3 P3, 

satisfy 

(10) T&T^p^fr, z = l,2,3. 

These functions define new potential f related to r by J5J). The new 
f-function satisfies dKP equation and due to condition (fTT)|) it fulfils 
condition (JEJ). As a consequence f satisfies equation 0. □ 

3.3. Explicit formulas for genus g = 0. In the case of projective 
line the function h reads 

_ (t-A 1 )(t-A 2 ) 
[ ) ~ (t- A 3 ) ■ 

Then, 

{M - A 2 ) (A 2 - A,) 

'^lA—Ajy 62 = Air 5 3=(A 3 -A 1 )(A 3 -A 2 ). 
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The formula © and the determinant formula (jBJ) allow to find pure 
7V-soliton solutions of the discrete ID Toda equation over finite fields. 

4. Examples 

Example 1. The vacuum solution of the discrete ID Toda in F 5 

First, in Figure ^ we present 0-soliton f -function taking values in 
the finite field F5. We fix parameters of the solution from F5 as follows: 
Ai = 0, A 2 = 2, A 3 = 3. For this example the Remark at the end of 
section |2~2"1 and formula Q imply relations Tff = (5 1 ) _2 f and T^f = 
(S 3 ) 2 f. Because 61 = 4 and 5 2 = 3 the period in n x variable is equal to 
4, and in n 3 variable is 8. It is also evident from the picture that the 
reduction constraint T 2 f = ^T{ l T 3 f where 7 = 2 is satisfied. Notice 
also there are no values in a vacuum solution. 

Elements of F 5 are represented by: 

■ -0, -1, -2, M-3, M-A. 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 



Figure 1 . The vacuum solution f(ni, n 2 , r^) of the dis- 
crete ID Toda equation for n 2 = and 1. Axes: rii di- 
rected to the right (from to 14), n 3 directed upward 
(from to 14). 

Example 2. The breather solution of the discrete ID Toda in F5 

The next example, presented in the Figure |2J is an analog of a 
breather solution taking values in the finite field F5. Parameters of the 
solution are chosen in a "symmetric" way from F 5 2 = F5 [x] / (x 2 + x + 1) 
i.e. a quadratic algebraic extension of the field F 5 by the polynomial 
(x 2 + x+ 1). An element ax + b G F 5 2 is denoted by (ab). The adequate 
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Frobenius automorphism is denoted by a F . Parameters are fixed as 
follows: 

A x = (00), A 2 = (02), A 3 = (03); C x = (10), C 2 = <t f (C x ) = (44); 
D x = (12), Ex = a F (Dx) = (41), D 2 = (24), E 2 = a F (D 2 ) = (32); 
The r-function is normalised by f(rix = l,n 2 = 0,n 3 = 1) = 1. Re- 
spective periods in variables n x and n% are equal to 12 and 24. The 
reduction constraint T 2 f = 7T{~ Tgf is satisfied (for 7 = 2) as visible 
in the figure. 

Elements of F 5 are represented like in the previous example by: 
■ - (00), - (01), ■ - (02), ■ - (03), ■ - (04). 




5 10 15 20 25 30 5 10 15 20 25 30 



FIGURE 2. The breather solution f(ni,n 2 ,n 3 ) of the 
discrete ID Toda equation for n 2 = and 1. Axes: n x 
directed to the right (from to 29), n 3 directed upward 
(from to 29). 
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